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Abstract 

We derive finite temperature versions of integral formulae for the two-point 
correlation functions in the antiferromagnetic XXZ chain. The derivation is 
based on the summation of density matrix elements characterizing a finite 
chain segment of length m. On this occasion we also supply a proof of 
the basic integral formula for the density matrix presented in an earlier 
publication. 
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1 Introduction 

Although many Yang-Baxter integrable quantum systems [11,31] with important ap- 
plications in diverse areas of physics have been discovered and intensely studied since 
the 60s, nowadays knowledge about the correlation functions of local observables in 
these models still mostly relies on approximate methods like the application of confor- 
mal field theories for critical models or direct numerical studies of finite systems. Any 
progress in the exact calculation of correlation functions is hindered by the complicated 
structure of the many-body wave functions, depending on parameters, the Bethe roots, 
which, in most cases, are not known explicitly. But also in those cases, where explicit 
expressions for the Bethe roots are available and the full finite temperature dynami- 
cal correlation functions could be calculated, namely for the impenetrable Bose [31] 
and Fermi [12, 13] gases, rather non-elementary mathematics is involved, and even the 
asymptotic form of the correlation functions is rather complicated. 

So far the only Yang-Baxter integrable system with a non-trivial distribution of 
Bethe roots in the ground state for which some more or less explicit exact results on the 
correlation functions could be obtained is the anisotropic Heisenberg antiferromagnet 
with Hamiltonian 

It acts on the tensor product of L spins |. The local operators a a , a = x, y,z, in (1) 
are the Pauli matrices and J > and A > — 1 are two real parameters, called exchange 
interaction and exchange anisotropy. 

Important progress in the calculation of the correlation functions of the XXZ chain 
was achieved in 1992, when Jimbo et al. [19] shifted the focus of attention from the 
two-point functions to the density matrix of a segment of length m of the infinite chain. 
From a physical perspective the density matrix of a sub-system seems a rather compli- 
cated object. Yet, it can be used to express the expectation value of any operator acting 
on sites 1 to m and, for Yang-Baxter integrable systems, has an interesting mathemat- 
ical structure, when generalized to include inhomogeneities on columns 1 to m of the 
six- vertex model associated with the Hamiltonian (1). As was shown later [1,2], the 
inhomogeneous density matrix elements satisfy a set of functional equations in the in- 
homogeneity parameters which was termed the reduced q-Knizhnik-Zamolodchikov 
equation (rqKZ equation). 

In [19,20] (multiple) integral formulae for the inhomogeneous density matrix ele- 
ments were obtained. These played a decisive role in the further developments of the 
following years. Their generalization to the finite temperature case proposed in [16] 
will be proven below and is the starting point for our derivation of integral formulae 
for two-point functions at finite temperature. Our reasoning will be based on the al- 
gebraic Bethe ansatz solution of the XXZ chain and will make extensive use of the 
Yang-Baxter algebra. Such type of approach was pioneered in [27], where the integral 
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formula for the density matrix of [19, 20] was rederived and a longitudinal magnetic 
field, not accessible by the original symmetry based approach of [19,20] was incorpo- 
rated into the calculation. 

At least at zero temperature the integral formulae for the density matrix turned out 
to be a means to derive completely explicit results for short-range correlation func- 
tions. This was first demonstrated in the isotropic XXX case [5] (A = 1 in (1)) with 
the example of a special matrix element of the density matrix, called the emptiness 
formation probability [32]. In succeeding articles [6,7,34] further explicit examples 
of short range correlation functions were obtained by calculating integrals, and the 
technique of [5] was extended to the anisotropic case [21,22,40]. The knowledge of a 
number of explicit examples then helped to develop an ansatz [8,9] for the general in- 
homogeneous density matrix element in the XXX case which did not involve integrals 
anymore. This was subsequently proved and further generalized and simplified [1,2,4]. 
Finally it was brought into a form which suits for taking the homogeneous limit [3]. 

The 'Bethe ansatz approach' of [27] was further developed in [24], where integral 
formulae for two point functions and for a certain one-parameter generating function 
of the S z -S z correlation functions were obtained. The important feature of these for- 
mulae is that they are different from those obtained by summing the integrals for the 
density matrix elements naively, in that they involve only a sum over as many (multi- 
ple) integrals as lattice sites, while the number of terms in the naive sum is growing 
exponentially. In a sense the formulae of [24] give an implicit solution of the problem 
of re-summing the naive sum into larger blocks. This probably inspired the work [26] 
on the explicit re- summation of sums over density matrix elements. Note, however, 
that the formulae of [26] differ from those in [24]. They seem to be simpler and led to 
new explicit results for A = i in [25]. 

This work is on the generalization of [27] and [26] to finite temperatures. It is 
based on previous joint work of two of the authors with A. Kliimper, where the quan- 
tum transfer matrix approach [28,29,38,39] to the thermodynamics of integrable lattice 
models was combined with certain combinatorial results [24,31,37] from the algebraic 
Bethe ansatz. As a result integral formulae for two-point functions [14] and density 
matrix elements [15, 16] at finite temperature and finite longitudinal magnetic field 
were obtained. In this work we shall present a proof of the integral formula [16] in the 
general case. We further apply the summation formulae of [26] to the finite tempera- 
ture density matrix and obtain new integral representations for two-point functions at 
finite temperatures. 

In the next section we review the algebraic Bethe ansatz approach to finite temper- 
ature correlation functions [14]. Section 3 is devoted to the density matrix. We explain 
its origin in physics and show how it naturally fits into the formalism of section 2. The 
remaining part of the section contains a proof of the multiple integral formula for the 
density matrix elements proposed in [16] and shown in theorem 1. Some technical 
parts of the proof are deferred to appendix A. In section 4 and appendix B we show 
how to sum up the density matrix elements into efficient formulae for the generating 
function of the S z -S z correlation function and for the correlation function (o^o m+1 )r./,. 
We conclude with a summary and an outlook in section 5. 
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2 Algebraic Bethe ansatz approach to finite tempera- 
ture correlation functions 

Before coming to the definition of the density matrix and to its representations as ex- 
pectation value of non-local operators in a certain Bethe ansatz state and as a multiple 
integral we have to fix some basic notation. We follow our earlier paper [14], where 
more details and proofs of most of the statements cited below can be found. 

In the framework of Yang-Baxter integrability the XXZ chain with Hamiltonian (1) 
is the fundamental model [11] associated with the solution /?(A,,/i) € End(C 2 ® C 2 ) of 
the Yang-Baxter equation 



R^MRa^MR^M = # py {ft, V )^«Y v )^a"p" A 1 ) (2) 



defined by 



(\ 0\ 

o HM c(M o 

c(X,/i) 

\0 1/ 



(3) 



KM = 



sh(A,— //) 



c(X,//) = 



sh(ri) 



(4) 



sh(X-Ai + Tl) sh(X-/i + Ti) 

For the exact calculation of thermal averages at temperature T we need to express 
the statistical operator of the canonical ensemble 

-Hxxz/T 



PL 



(5) 



in terms of the 7?-matrix (3) in a way that allows us to utilize the Yang-Baxter algebra 
(see (13) below) efficiently. For this purpose we introduce an auxiliary vertex model 
defined by a monodromy matrix 

W = (cm om) . = ^>I)*M-^) ..^fiKHU) • ( 6 ) 



VC(X) D(X)) j 

This monodromy matrix acts non-trivially in the tensor product of an auxiliary space 
7 = 1,.. . ,L and € 2N 'quantum spaces' 1, . . . ,N which are all isomorphic to C 2 . By 
t\ we mean transposition with respect to the first of the two spaces into which R(k,ij) 
is embedded. The parameter |3 is inversely proportional to the temperature T. By 
definition, 

27sh(ri) 



P 



(7) 



The monodromy matrix (6) has many remarkable properties. First of all it can be 
used to express the statistical operator (5) as a limit, 



p L = limtrj #7i(0)...7i(0), 



(8) 
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with the implicit identification A = ch(r|). The number N in (8) is called the Trotter 
number, and the limit is referred to as the Trotter limit. 
For large Trotter numbers the operator 

p^ = tr L ^7i(0)...r L (0) (9) 

is a good approximation to the statistical operator. The thermal expectation value of a 
product of local operators acting on sites j,...,j + m—l of the spin chain is then 
approximated by 

, (1) (m) v tri..xpjy,LAj ■■■* j+m - l 

\ j ■ ■ -Xj+m-lJNJ = 

J J ui...LpN,L 

_ trj J v(trr(0)) J " 1 trX( 1 )r(0)...trXHr( )(trr(0)) L ^ J ' +1 

tr L .^(trr(0)) L 

Here the expression 

Z n ,l = tri...L Pn,l = tr T ^ (trr(0)) L (11) 

in the denominator is the finite Trotter number approximant to the partition function 
Zi = tri. z,exp(— fl^/r) of the XXZ chain of length L. 

The operator occurring on the right hand side of equation (11), 

t(X)=trT(X)=A(X)+D(X), (12) 

is called the quantum transfer matrix. By construction it can be diagonalized by means 
of the algebraic Bethe ansatz. This follows [31] since the monodromy matrix (6) 
defines a representation of the Yang-Baxter algebra, 

R Jk (X,fi)Tj(X)T k {fi) = T k (iS)Tj{\)R jk {\^) , (13) 

and has a pseudo vacuum |0) = [(J) ® ( J )] ® NI2 on which it acts like 

C(X)\0)=0, 

A(X) |0> = m ( - | ,X) |0> = a(X) |0> , D(X) |0> = ^ (x, J) |0) = d(X) |0) 

Then (see [31]) the vector 

\{X}) = \{X j }f =l )=B(X l )...B(X M )\0) (15) 

is an eigenvector of the quantum transfer matrix t(X) if the Bethe roots Xj, j = 1 , . . . ,M, 
satisfy the system 

a(Xj) shfo-Att + TQ 

d(Xj) l^sHXj-Xk-r]) K } 

Hi 



(14) 



6 



of Bethe ansatz equations. The corresponding eigenvalue is 



Note that the 'creation operators' B(kj) mutually commute. Therefore the eigenvector 
(15) is symmetric in the Bethe roots, whence our notation | {k}}. 

Equations (16), (17) are a means to study the spectrum of the quantum transfer ma- 
trix (12). It is easy to see that for every finite Trotter number N and high enough tem- 
perature T the quantum transfer matrix has a non-degenerate eigenvalue with largest 
modulus at X = which is characterized by a unique set {X} = {Xj}^/^ of M = N/2 
Bethe roots that all go to zero as T goes to infinity. We call it the dominant eigenvalue 
and denote it by Ao(A,). The corresponding eigenvector \{k}) will be called the dom- 
inant eigenvector. It is a crucial property of the whole construction that the dominant 
eigenvalue alone determines the thermodynamic properties of the XXZ chain in the 
thermodynamic limit L — > oo. This can be seen as follows. The free energy per lattice 
site in the thermodynamic limit follows from (1 1) as 

f(T) = -T lim - lim lntrj N t L (0) . (18) 

But the two limits in (18) commute according to [38, 39], hence 

f(T) = -T limlnAo(O). (19) 

iV— >°° 

Applying a similar reasoning to the expectation value of the product of local operators 
(10) in the thermodynamic limit we obtain [14] 

(X W X W v _ lim ({M|trX( 1 )r(0)...trXWr(0)|{M) 

where \{k}) denotes the dominant eigenvector. This means that the dominant eigen- 
vector determines all thermal correlation functions of the system in the thermodynamic 
limit. In other words, it determines the state of thermal equilibrium completely. 

The expression under the limit on the right hand side of (20) expresses the finite 
Trotter number approximant of the expectation value of the product of operators on the 
left hand side of the equation as an expectation value of monodromy matrix elements 
in a Bethe state. But this is still not good enough for applying the Yang-Baxter alge- 
bra (13) since the arguments of the monodromy matrix elements are all zero, and the 
relations in the Yang-Baxter algebra for commuting two monodromy matrix elements 
with the same spectral parameter X are singular. Therefore we first regularize the ex- 
pectation value by introducing a set of mutually distinct 'inhomogenei ties' £i, . . . 
and then calculate the correlation functions by means of the formula 

j=l,...,m 
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Since the Hamiltonian (1) commutes with the z-component 

L 
7=1 

of the total spin, thermodynamics and finite temperature correlation functions can still 
be treated within the 'quantum transfer matrix approach' described above if the system 
is exposed to an external magnetic field h in z-direction [14]. The external field is 
properly taken into account by applying a twist to the monodromy matrix (6), 

T(l)^T(l)[^ T J /1T ). (23) 

In [28, 29] techniques were developed to treat the Trotter limit in (19) analytically. 
The main idea is to associate with every Trotter number /V G 2N a meromorphic 'aux- 
iliary function' 

</(*,) ffi shfl,- ^ + tQ 
a^(A) = — Y-M-— — — -, (24) 

where the Xj are the Bethe roots that characterize the dominant eigenvalue, and to 
consider the sequence of these functions in the complex plane, rather than sequences 
of sets of Bethe roots. This is recommendable, since the distribution of Bethe roots for 
large Trotter number N is not smooth but develops a limit point at X = 0. The sequence 
of functions <\ N , on the other hand, converges to a limit function o which is analytic far 
enough from X = 0. The strategy is thus to express all physical quantities of interest 
as contour integrals over on contours far away from the origin and then take the 
Trotter limit. 

In particular In a# can be expressed this way [14], which yields a non-linear integral 
equation for a#, 



h 

Inojv(X) = --+ln 



N 
1 



sh(A,-g)sh(A,+ ft+Ti) - 

sh(X+g)8h(X-|+Tl). 

do sh(2r|)ln(l + ajv(to)) 



-L 



c 2ni sh(A, — co + ri) sh(A,-(o-ri) 



(25) 



Notice that we included the magnetic field here according to the recipe explained 
above. The integration contour C depends on r\ (see figure 1). As is very impor- 
tant, however, it does not depend on N. The Trotter number /V merely appears as a 
parameter in the inhomogeneity of the integral equation. Hence, the Trotter limit can 
easily be performed. Using also (7) we obtain 

lnam _ h 27sh 2 (rQ r da sh(2rQln(l +a(co)) 

{) T rsh(X)sh(X + ri) J c 2%i sh(X- co + r|) sh(X- co-r|) ' 

A physical interpretation of the auxiliary function a stems from the 'free Fermion 
limit' A — > (see e.g. [17]), where the Hamiltonian (1) describes the highly degenerate 
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Figure 1 : The canonical contour C in the off-critical regime A > 1 (left panel) and in 
the critical regime |A| < 1 (right panel). For A > 1 the contour is a rectangle with sides 
at ±i| and ±y, where y is slightly smaller than ^. For |A| < 1 the contour surrounds 

the real axis at a distance |y| slightly less than 4^. 



XX chain. In this limit the function 1/(1 + a) becomes the Fermi function. Thus, 
we may think of 1/(1 + a) as generalizing the Fermi function to the interacting case. 
Thinking in terms of Fermions it is natural to set a = 1 / a, and interpret 1 / ( 1 + a) as 
the Fermi function for holes satisfying 1/(1 + a) + 1/ (1 + a) = 1. 

In the thermodynamic limit the free energy per lattice site as a function of temper- 
ature and magnetic field can be expressed as an integral, involving either the auxiliary 
function o or its reciprocal a, over the same contour as in (26) [14, 30] 

f ( hT )= h T [ d& sh(Ti)ln(l + q(<p)) = h r d& sh(rQln(l + a(a))) 

J{ ' ; 2 J c 2%i sh(a>)sh(a> + Ti) 2 J c 2%i sh(co) sh(ro-r|) ' 1 ; 



3 The density matrix 

The density matrix is a means to describe a sub-system as part of a larger system 
in thermodynamic equilibrium in terms of the degrees of freedom of the sub-system. 
Roughly speaking it is obtained by taking the (normalized) statistical operator and 
tracing out the 'unwanted' degrees of freedom. It can be defined rather generally, but 
we content ourselves here with the example when the sub-system consists of the first 
m sites of the XXZ chain of length L> m. Since we want to include a magnetic field 
in z-direction into our considerations we have to modify the statistical operator (5) as 

p L = exp(-(H xxz -hS z )/T). (28) 

Then the density matrix of the sub-system consisting of the first m lattice sites is de- 
fined as 

D L (T,h) = tVm+l - LpL . (29) 
tri...LpL 
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By construction, the thermal average of every operator A acting non-trivially only on 
sites 1 to m can now be written as 

\ A )T.h = - — = : =te\..mM...mPiXT,h), (30) 

tri..xpL tr L .. L pL 

where A\ w is the restriction of A to a chain consisting of sites 1 to m. In particular, 

every two-point function of local operators in the segment 1 to m of the XXZ chain 

can be brought into the above form. 

In order to make contact with the formalism explained in the previous section we 

shall refer to a particular local basis. For our purpose it is most convenient to choose 

the gl(2) standard basis consisting of the 2x2 matrices e?, a, (3 = 1,2, with a single 

non-zero entry at the intersection of row (3 and column a. The canonical embedding 

of these matrices into the space of operators on the space of states (C 2 )® L of the 

Hamiltonian (1) will be denoted by ef^, j = 1, . . . ,L. For their restriction to the first m 

sites we use the same symbols. Then we find the expressions 

D L ^(T 1 h) = tri.. M ei^...e m ^D L (T 1 h) = . . .e m ^) T , h (31) 

for the matrix elements of the density matrix, which hold for all L and are still valid in 
the thermodynamic limit. To denote the density matrix in the thermodynamic limit we 
simply leave out the subscript L in (3 1). 

For L — > °° the right hand side of (31) is of the form of the left hand side of (21) 
with 7 =1 and X^> = et 1 . Hence, it can be calculated as 

«' ■ -'«Z)t,h = h ,..^ DW ££fo- ' ' ' «-> ' <32) 



where 



w a,.. a _ ({MI^(^)-r- :( UIW) 

p ,.. Pm Ki,...,w • <{x}|ny =1 /(^)im> 



(33) 



is the 'inhomogeneous finite Trotter number approximant' to the density matrix el- 
ement (31) in the thermodynamic limit. In order to simplify the notation we have 
suppressed the dependence on T and h on the left hand side. Notice the formal simi- 
larity between (29) and (33). In both cases the denominator is the trace of the term in 
the numerator. 

In [16] a basic integral formula for the inhomogeneous finite Trotter number ap- 
proximant (33) of the density matrix elements was proposed. The formula involves 
certain positive integers aj and |3j derived from the sequences of indices (cc n )™ =1 , 
($n)n=i characterizing the matrix element. The indices take values 1,2 corresponding 
to spin-up or spin-down. We shall denote the position of the jth up-spin in (oc M )™ =1 
by aj~ and that of the Mi down-spin in ((3 M )^ =1 by (3^. Then, by definition, a+ = 
ccj^+|_^. +1 , where |a + | is the number of up-spins in (a„)™ =1 , and j = 1, . . . , |a + |. Sim- 
ilarly (3j = P^_| a+ i, k = \a + \ + 1, . . . , |a + | + ||3~|, with |(3~| denoting the number of 
down-spins in (Pn)^Li- 
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Theorem 1. Due to the conservation ofS z the matrix elements D^pJ.'."^ 1 (^1 , • • • , \rn) 
vanish ?/|cx + | + |P~| 7^ m. For |a + | + |f$~| = m they are given by the formula 



fl <K) S::.tK..-,w 



n 



P7- 1 



j i! / ,..„,, . — r , \ EI sh(co 7 --^ + Ti) n sh(co y --^) 
. J =|a+|+i y c2ra(l + a^(a) J )) £=i fc= |i +1 

det(-G((o 7 ,^)) 



ril< 7 <fe<m sh (^ - y Sh(t0 ; - -C0fc-T|) 



(34) 



The function G((fl, ^) depends on the Trotter number through o-n- It has to be calculated 
from the linear integral equation 



G(U) = 



sh(ri) 



sh($-X) sh(^-X + ri) 



+ 



/ 2ro(l 



J(0 



8h(2Tl)G(Q)^) 



/ c 27ci(l + ajv(a))) sh(X-a> + r|) sh(X-a>-r|) 
77ze contour C in (34), (35) is the same as in figure I. 



(35) 



Note that the function G((0, introduced in [14] generalizes the 'density function' 
p(A,), that determines the ground state energy and the zero temperature magnetization, 
to finite temperatures and to the inhomogeneous case. Taking the Trotter limit and the 
homogeneous limit (for the latter compare [27]) in (34) we obtain an integral formula 
for the density matrix elements (32). 

Corollary 1. Integral formula for the density matrix of a segment of length m of the 
XXZ chain in the thermodynamic limit. 



(eifa---emf£) Tjh - EI 

m 

n 



d(Dj 



j=\Jc 27ri(l + a(Cfl;)) 

d(tij 



.. =| ^ +1 Jc27ri(l + a(t0;)) 

af- 1) G(to y ^)|^ =0 



det 



sh^'- 1 ((o J --ri)sh ra - 5 i((o J 
shP7- 1 ( (0j . + n)sh ra ^7((o ; ) 

1 



(k-iy. \Hi<j<k<mM®j-®k-n) 

Here G((0, is defined by (35) with <xn replaced by a. 



(36) 
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In the following we shall give a proof of theorem 1 . The proof is based on combin- 
ing the two following lemmata. Lemma 1 is new and of its own interest. It describes in 
a compact form the general left action of a string of monodromy matrix elements T^fc) 
conserving the z-component of the total spin on a state of the form (0|C(A,i) . . .C(Xm)- 
Lemma 2 is a simple corollary of a result originally obtained in [14]. 

Lemma 1. The general left action. For m,M e N choose Xi,. . .,XM+m £ C mutually 
distinct. Define 



a(X) sh(X-X k -r\) 
and for j = 1 , . . . , m and £j G { 1 , . . . , M + m) 

|£ + |=card{l,...,M}n{^ =1 

and = Xm+]- Then 
(0\C(X l )...C(X M )T^(^)...T^(U 



(37) 



(38) 



M+m 



I_ (-sh(Tl)) 

£\ ,...,£ m — 1 
tj¥*kfi>rj¥* 



\t+\ 



M 



I] b M (X £j ) Y\a(X ej )l] 



1 



■j=\a+\+\ 



n b (h v xtj) 

\<j<k<m 



1 



n n sh(^.-y 

k+M^ij 



n n sh^.-^-T!) n m^-w 

;=1 fe=l 



P7 



n n sh (H -^+ T i) n m^-w 

■j=\a+\+l k=l k=$j + l 



(0| 



M+m 

n 

£=1 



(39) 



A proof of this lemma is given in appendix A. Note that the Bethe ansatz equations 
are not used in the derivation of (39). However, (39) also applies if (0|C(A,i) . ..C(Xm) 
is the dominant (left) eigenvector of the quantum transfer matrix. Then necessarily 
M = N/2 and bM = cijv- This is the case needed in the proof of theorem 1. 

After inserting (39) into the right hand side of (33) it reduces to a sum over matrix 
elements of the form 

(ft + }u{x-}|{x»=(oi[ n c & + )}\ n c(A,-)i[n^)li o >> w 

where {^+} c {^}, {X~} c {X} and card{^+} U {X~} = card{X} = N/2, i.e. we have 
to deal with partitions of sets {X} into disjoint subsets, for which we shall adopt the 
general notation p n {X} with n denoting the number of subsets. We further frequently 
use the shorthand notation |X| = cardX for the number of elements in a set {X}. 
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Lemma 2. The matrix element formula. Let {X} be the set of Bethe roots that para- 
meterize the dominant eigenvalue of the quantum transfer matrix, and let {^} be a set 
of mutually distinct complex numbers from inside the contour C such that {^} H {%} = 
0. Then, for ({£+},{£-}) G p 2 {^} and ({X+}, {X~}) G p 2 {k) with = 



<{^+}U{A,-}|{A,» 

<wim>n?=iAo(^) 

n 



| = V(^)(l + a^T)) 



}J a(X+)a^(X+) 



detG(A+,# 



det 



(41) 



where M = N/2 and where the prime denotes the derivative with respect to the argu- 
ment. 



This lemma follows by applying elementary manipulations to the formula exposed 
in lemma 1 of [16]. 

Let us now proceed with the proof of theorem 1 . In order to be able to combine 
(39) and (41) into the right hand side of (33) we must distinguish between Bethe roots 
and inhomogeneities in the sum in (39), i.e. we have to split it into contributions with 
summation indices running from 1 to M = N/2 and from M+ 1 to M + m, respectively. 
Then 

M M+m 

£ £ £ (42) 

({E+},{e-})€ P2 (z m ) £ E p...,e E +_ n =i e e -,..j E -=M+i 

Here Z m = { 1 , . . . , m} and n = card{£~ } by definition. We shall employ the convention 
to enumerate the elements in {e^ 1 } in increasing order, 

£+<£+<•••< e+_ n , Ei < £ 2 <•••<£„ . (43) 

Taking into account that £ + G { 1 , . . . , M} and £ - G {M + 1 , . . . , M + m} we see that 

j j 

we must set 

& + } = {\}7=* > {\-}={\.}u- ^ 

We enumerate the elements in {k + } and {^~} in such a way that 

The enumeration of the elements in {^ + } and will be fixed later, when it is 

needed. Combining now (33), (39) and (41) and using (42), (45) we obtain after some 



M+m 

E 

lj+l k for j+k 
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trivial cancellations 



M 

E 



({£-},{e+})ep 2 (Z m ) e+,...,e+ =1 £ =M+1 

t l *m-n 1 « 



M+m 

£ det[-G(V + ,g-)] 



7 fc 





e+<|cc+| 7 



n 



n F Ej(\-) 

n - - L — 



n - 

(46) 



Ej >|or 



eJ<|or 



\<j<k<m 



mm j 

n n sh(^.-y 



det 



sh(A.,; 



e ; . >|a 



Here we have used the Bethe ansatz equations in the form Oat(A,£ + ) = — 1, and we have 
introduced the shorthand notations 



ocj-i 



Fj(a>)= n sh(to-^-Ti) f[ sh(©-y, j = l,...,|a+|, (47a) 



~ -i 



k=aj+l 



~j m 

F j(®) = EI sh(co-^t + Ti) f] sh («-^), j= |a+| + l,...,m. (47b) 
The last two terms in (46) can be further simplified. We first rewrite them as 



1 



n n sh(^-y 

k+M^ij 



det 



l 



i 



l<Mm-n Sh (\ t -\ f ) 



k 
1 



1 



\<j<k<m—n 



m—n n 
n m—n ^ 



n - 

l<j<k<n M\-~h £T ) 



nn sh( ^_^ 



n _ 

\<j<k<n 



(48) 



Then two permutations P,Q e 6 m exist, such that 



£, = ^i\ J = 1, 



eJ = P(n + j), j=l,...,m-n, (49a) 
^J" = ^fi(n+y)' 7=l,---,w-n. (49b) 
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After inserting (49) into the right hand side of (48) and shifting the multiplication 
indices the products can be combined into 



1 



n n sh(^-y 



1 



det 



sh(^ 



•k n 



U 

l<j<k<m 



1 



sh(A,< 



pk - ^tpj 



ish(£ 



Qj 



sign(PQ) 11 — -- 

l</i<m Sh ^4-^) sh (^-^) 



£>Qk) 

(50) 



Inserting this back into (46) we conclude that 



({e+},{e-})e P2 ( 



i-n F e+ fte 



FT -L 2 



\<j<k<m 
M 

E 



M+m 

E 



sign(PQ)det[-G(^ + ,y-)] 
2 

i Ui<j<k<mMh.j-h k -n) 



t -& - for j^k 

j k 



+ a N (X f 



n i+oa 



+ oat(A,^ _) 



ej<|cx+ 



ej>|or 



eJ<|or 



£j >|or 



(51) 



Note that we could leave out the exclusions in the sums over £ £ + , j 



i,...,m — n, 

since det[— G{Xi vanishes if two of the £ £ + coincide. Note further that Q is not 

£ j j 
entirely defined by (49), since so far we did not fix a convention for enumerating the 

£+. It follows from (45) and (49b) that 



(52) 



(53) 



^ = Ht =Z=Qj = \--M & Qj = £ e j-M, j=l,...,n. 
In order to fix Q uniquely we stipulate that 

Q(n+1) < Q(n + 2) < ■ ■< Qm. 

This then also fixes the enumeration of elements in {^ + } through (49b). Our freedom 
in the choice of the enumeration is explained by the fact that sign(<2) is multiplied by 
det[— G(k,£ + ,^|~)]. The product is totally symmetric in the 

In the next step we would like to make the dependence of the terms under the sum 
in (51) on the inhomogeneities more explicit and rewrite the sum over the t z - in a more 
convenient form. The inhomogeneities occur only in the last two products, where we 
can directly use (52), and in the product in the denominator, which we shall write as 



i<j<k<m \<j<k<m 



CO 



,00 



(54) 



£j -M 
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In the multiple sum over £ £ - , . . . , £ £ - all summation variables must take mutu- 
ally distinct values. For this reason each term in the multiple sum defines a partition 

({5+},{5-})ep 2 (Z m ) via 

{5-} = {£ e --M,...,£ -M}, {5+} = Z m \{5-}. (55) 

1 n 

By our convention the elements in {S 3 * 1 } are ordered, 8^ < 8^ if j < k. Hence, a unique 
Re & n exists, such that 

8 Rj = £ £ ,-M, 7 = 1,. ..,n. (56) 

On the other hand, each pair consisting of a partition ({8 + },{8~}) £ p2(1^ m ) with 
|8~| = n and a permutation R £ ©" uniquely determines an n-tupel -, . . . ,£ e -) of 

1 w 

summation variables by equation (56). We conclude that with the substitution (56) the 
sum over the £ £ - in (51) may be split as 



j 



M+m 

£ ... = £ £ •••• (57) 

V'-'V =M+1 ({8 + },{5-})e P2 (Z m ) *e6" 

£ _^ _ for I 8 l=" 

The permutation R E & n has a natural embedding into 6 m defined by 

R=[ Rj if J =1 — n (58) 
(j i f j = n+l,...,m. 

We would like to split off 7? from the permutation Q defined in (49b), (53). Using (52), 
(56) and (58) we obtain 

Qj = £ £j -M = 8~. & 8j = QR l j, 7=1,. ..,n. (59) 

It follows that Q = QR 1 satisfies 

Ql <Q2< - - <Qn, Q(n + 1) < Q(n + 2) < ■ ■ ■ < Qm. (60) 

Thus, there is a one-to-one correspondence between the permutation Q and the parti- 
tion ({8+}, {8-}) £ P2 (Z m ), 

87 = 67, 7=1,. ..,n, 8t = 2(71 + 7), 7 = l,...,m-7i, (61) 
and 2 does not depend on 7?. With (49b) we find that 
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Now we are in a position to rewrite the sum (5 1) in a form suitable for transforming 
it into a multiple integral. We use (57) in (51) and insert (52), (54), (56) and (62), 

\<j<k<m 

M _ 

= 1 £ I sign(P<2) 

({£+},{£-})e K (Z ffl ) U,...,€ e+ =1 ({8+},{8-})e K (Z ffl ) 
1 m ~" |8-|=n 

sign(i?)det[-G(^ ^ 5 +)] 



*e6» Ul<j<k<m sh((Oj — CO* — Tl) 



m-n 

n 



;=1 a^(^ e+ ) 

£+<|a+| 1 



m-n F z+ V^L+ 

n - -j- 



y 1+0, 



=1 l +Mt>8-) 



ej>|cT 



Zj >|cr 



(63) 



Recall that we have assumed that the inhomogeneities £,j are mutually distinct and 
distinct from the Bethe roots. Hence, our canonical contour C can be decomposed into 
two simple contours <B and I , such that 



C = $ +1 



(64) 



and <S contains only Bethe roots and no inhomogeneities, while / contains only in- 
homogeneities but no Bethe roots. Then 1/(1 + a#(a))) is holomorphic inside / and 
meromorphic inside <B with only simple poles, located at the Bethe roots [14]. On 
the other hand, — G(cfl,£y) considered as a function of (0 is holomorphic inside <B and 
meromorphic inside I with a single simple pole with residue 1 at t,j [14]. Thus, 

\<j<k<m 



M 



III sign(PQ) 

({£+},{£-})e K (Z m ) =1 ({8 + },{8-})ep 2 (Z m ) 



n 

■7=1 



d(0 o 



I sign(J?)[-G(a> er ,^- )] ... [-Gi^^-J 



i 2%i 

det[-G(^,^ 5 +)] 



s 



det[-G(cO E -,^ 5 - 



[] sll(CD 7 - — COjfc — Tl) 
y l<j<k<m 



-n F e+(h 



n 7 



y =1 a^(^ g+ ) 



m-n ~F £ + ( ) 

n - ^ 



}J l + Ma> eT ) 



£j<|or 



ej >|cr 



f = \ l + ayv(w e7 ) 



£■ >|cr 
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I 

({£+},{£-})ep 2 (Zm) 



m— n p d(fl £ + 1 I" n p d(fl £ 



Ill<y<Jt<mSh(Q) y --Q)t-Tl) 



dtflj 
C 2jff 



({8+},{5-})e P2 (Z ffl ) 
|8"|=n 



sign(PQ) det[-G(to e7 , )] det[-G(to e| , £ 5 + )] 

J K J K 



i—n _F g +(CO e +) 



JJ 1 + Mro e+ ) 



=det[-G{<Oj&)] 

m—n F c + f(0o+ ) 

n i j 

;=1 l + M(D e +) 

e+>|a+| 



L }i l+M(D e7 ) 



!J l + ayv(ro eT ) 



7=1 



eT<|a+| 



>|or 



| = Vc 27ii 1 + ajv (<*>;) 



n 



det[-G(oo ; ,^)] 
Y[ sh(a>_/ — to*: — n) ' 

\<j<k<m 

(65) 



and our proof of theorem 1 is complete. Note that we used the Laplace expansion for 
determinants in the third equation. 



4 The two-point functions 

In this section we shall derive a number of formulae for two-point correlation func- 
tions. Since the derivations concern only the 'algebraic part' of the integrals, namely 
the product 1/ III </<*</» sh((Oj — Cfl£ — T|) and the functions Fj, Fj introduced in equa- 
tions (47), we may strongly refer to the zero temperature results obtained in [26]. 

Let us start with a certain one-parameter generating function of the S z -S z correla- 
tion functions that was introduced in [18]. It is defined as the thermal average 

m 

4>(cp|m) = <exp{(p £ ej\}) T h . (66) 

7=1 

The function 4>(cp|m) generates the two-point function through the action of a second 
order differential-difference operator, 



(S\S z m ) Th = l(2D 2 Jl-4D m d (? + l)*(q>|«) 



(67) 

9=0 



Here D m denotes the 'lattice derivative' defined on any sequence (z n )n€N of complex 
numbers by D m z m = z m - z m -i- 
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The function 4>((p|m) is particularly useful and comparatively simple, since it is 
closely connected to a twisted version of the quantum transfer matrix as can be seen 
by inserting = (q e f ) for n = 1, . . . ,m into our basic formula (21). Then 

*(<p|m) = lim lira *W(<p|$i, • • • ,U , (68) 

where 

*<»>(»,£, E ) . <WIIg..(Aft/)+^))IW) 

(W | W )n™ =1 A„(y ' <69) 

is a polynomial of degree m in e^. 

We see from (34) that the Trotter limit and the homogeneous limit £,j — > com- 
mute. Hence, it makes sense to perform the Trotter limit in (34) and to consider 
'inhomogeneous generalizations of correlation functions'. An example for such an 
inhomogeneous generalization is the density matrix element itself, 

Df i:: ^ u . . . ,u = imDwj;:;^,. . . ,u , (7 o) 

obtained from (34) by taking the Trotter limit, i.e., by replacing a^v by a. 
Let us define 

f n (z >l ,...,U= I 0g::a&,-,$»). C 71 ) 

ai,...,a m =l,2 
h0t m =2ra— n 

Then 

m 

4>(cp|^,...,U= lim$W((p|^ 1 ,...,U = £e( m -") ( PA(^ 1 ,...,U (72) 



«=o 

is the inhomogeneous generalization of 3»((p|m), equation (66). Following [26] we 
shall derive an integral representation for the coefficients / n (£i, . . . , £ m ) of the polyno- 
mial 4>((p|^i , • • • , ^ m ) • This gives an integral formula for the generating function (66) in 
the homogeneous limit. The coefficients f n (in the homogeneous limit) have a simple 
physical interpretation. They represent the probability for a segment of length m of the 
infinite chain to be in a state with z-component of the total spin equal to n — m/2. For 
instance, for n = m we find f m = lim<p_ s ._ 00 < I>((p|m) which is the probability to have 
maximal S z = m/2 and which is, of course, equal to the emptiness formation proba- 
bility for spin-up. Thus, one may think of the f n as generalizations of the emptiness 
formation probability. 



Theorem 2. Define 



= -o J^rr e \ (73) 

v ' sh(co-^)sh(co-^ + r|) 
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and for n = 0, . . . , m 



tfa,aij) for\<j<n, 
t((Oj,Z,k) forn<j<m. 



(74) 



Then the following integral formula holds for the coefficient f n in the expansion of 
the inhomogeneous generalization of the generating function of the S z -S z correlation 
functions with respect to e*, 



/n(^l) • • • j^m) — 



n\{m — ri)\ 



n 

■y=i 



d&j 



c 2jti(l + a(cfl y -)) 



n 

j=n+\ 



d(tij 



c 2ra(l + a(ro J -))_ 



fj fr sh K -$k-n) sh((Dj - jU) p. sh^-^+rpsh^-^) 

det(-G(cfl;,^)) det(M M (co;,^)) 



ni<;</t<m M^k ~ Ul<j<k<m M^j - kk) 

A proof of this theorem is shown in appendix B. 
Corollary 2. Integral formula for f n in the homogeneous limit. 



(75) 



fn 



" r d(Ojsh m ((Dj-r[)sh m ((Dj) 

}J Jc 27ri 1 + 0(0);) 

(_1 yj+m(m-l)/2 



n\{m — ri)\ 



' ™ y J(o ; sh m (ro ; +ri) sh ra ((0 J ) 
Jc 2ri l + a(a>;) 

i m i 

n 



JJ JJ sh(CQ; -CO fc -Tl) yi^j sh(CQ; - & k + T\) 



det 



det 



(76) 



(*-l)! 

Similar but slightly more complicated formulae can be derived for the general two 
point functions for points at distance m defined by 



G®l(m;T,h) = lim G»g(£i, . . • ,£m+l) , 



where 



m—l 



Gp§(^l 5 ---^m+l) — J! ^np5(^l,---,^m+l), 



<?np§(^l7 • • • )^m+l) — 



a2,...,a„,=l,2 
a 2 H ha m =2m-n-2 



r na 2 ...a m y ( p e \ 
^pa 2 ...a m 5^1'---'WlJ- 



We shall consider 
as an example. 



(o+o- +1 ) T , h = Gii(m;T,h) 



21 , 



(77) 

(78) 
(79) 

(80) 



20 



Theorem 3. The coefficients grv^i^x-, • • • , ^m+l) in the expansion (78) of the inhomoge- 
neous generalization of the two-point function (o^o~ +1 )T : h are given by the integrals 



gnnfciT • • >£m+i) — 
(-1)" 



n\(m — n — \ )\ 



n+1 

n 



J(0, 



m+1 

n 



•=i A: 2ra(l + a(co 7 -)) J L^Vc 27li (l + a K)) 



J(0, 



det m+ 1 ( G(co, , y ) det ra _ i (M n ( to y+ 1 , ^+ 1 ) ) 



Ill<;<fc<m+lSh(^-^) IL<7<Jt<mSh(^t-^) sh(tO m+ i - tOi + T|) 

nr=inj!sh(a) J -^-Ti)n7i +2 sh(a) 7 -^+Ti) [*h » 

nri sh K-^) 



UT=2 U%\ M®j -COjfe-Tl) n7+; +2 sh(C0; - (0* + T]) 



m+1 



£=2 y=2 



(81) 



A proof of the formula is easily obtained by combining the ideas of appendix B 
with those of [26]. 



Corollary 3. Integral formula for in the homogeneous limit. 



„ 21 

g«12 



n+1 , 
.j=\JC 



'"+' ' Jro J sh m ((o J --ri)sh' n (a) J -) 
c 27ti 1 + o(t0 7 ) 



m+1 



n X 

j=n+2 Jc 



J(o J sh ra ((o J +ri)sh m (a) ; ) 



2ra 



( i \n+l+m(m— 1)/2 I" m n+1 

mi 



n\[m — n — \ )\ 



1 



m+1 

n 



l + a(a> y -) 
1 



det m+ i 



^ A 2j A =i sh(ro J - — CO* — ti) .4^2 sh(ro J - -co fc + Ti) 
det m _i 



a^-^G(to y ^)|^ =0 



(Jfc — 1 



(*-l)! 
1 



sh m (a)i) sh ra ((o m+ i) sh(co m+ i -coi + r|) 



(82) 



5 Conclusions 

We presented a proof of the integral representation of the density matrix (34) and 
obtained efficient formulae for a one-parameter generating function of the S z -S z corre- 
lation function and for (o^o m+l )r,h- The latter formulae were obtained by applying 
the re-summation technique developed in [26] in the finite temperature case. This 
technique applies to arbitrary finite-temperature two-point functions. 

We hope that our results will prove to be as useful as the multiple integral formu- 
lae in the zero temperature case and will eventually lead to explicit expressions for 
short-range correlation functions reaching beyond our preliminary numerical calcu- 
lations of the integrals in [10] (for the most recent progress at T = see [35, 36]). 
We believe that the integral formulae can be efficiently used to obtain high-order high- 
temperature expansions for the two-point functions, closely following the example [41] 
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of the emptiness formation probability. We further hope that we may be able to extract 
the long-distance asymptotic behaviour of the correlation functions directly from the 
integrals (for the zero temperature case compare [23, 33]). These and many other in- 
teresting questions are still open, but we may now be close to their answers. 

Acknowledgement. The authors would like to thank H. E. Boos and A. Kliimper for 
helpful discussions. This work was partially supported by the Deutsche Forschungs- 
gemeinschaft under grant number Go 825/4-2. 
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Appendix A: The general left action 



Here we shall provide a proof of lemma 1 . This means that we have to calculate the 
left action of a product of monodromy matrix elements on a state of the form 



M 



<{W=ii=<°in c (**) 



(A.l) 



where A4, . . . ,A,m G C are mutually distinct. We shall restrict ourselves to products 



7p m (^ m ) which contain as many 5s as Cs, i.e. products for which 



(A.2) 



Only these products have non-zero expectation value ({A}]^ 1 (£1) . . .T^ m (^, m )\{X}}. 

In order to obtain the action of a whole product of monodromy matrix elements 
we first of all need compact expressions for the action of a single monodromy matrix 
element. As was observed in [27] the known expressions for the single action (see 
e.g. [31]) can be written more compactly if one introduces the notation £y = Xm+j for 
the inhomogeneity parameters. This idea was extended to its full beauty in [16], where 
the equations 



1 



M M+l rM+l 

(01 n c (^) m^m+i)= e«(ww,m n , n 



M+l 

(0| HC(l k ), (A.3) 

k=\ 

k^l 



M M+l rM+l 

(01 ijc(i k )\D(i M+l ) = e m^+i) n 



1 



M+l 

1 n c (^)' ( a - 4 ) 
k=\ 

k^i 



M M+l M+l [M+l 1 

(01 i\c(i k ) b(w)= E E ^M^w) n 7^ — 

l k=l J £ l = l£ 2 =l Ljfc=l "Mi^J 

£2^1 



1 



■ M+l 

n fc^,^) 



M+l 

(0| n C(A*). (A.5) 



were obtained. 

These equations are an appropriate starting point for our proof. The natural thing 
to do is using them iteratively. Then each operator A or D generates a single sum and 
each operator B a double sum. The number n of these sums after acting with the first 
j monodromy matrix elements . . . T^ J is equal to the number of up-spins in the 

subsequence (cc n )^ =1 plus the number of down-spins in the subsequence {%) J n=l . This 



j 


1 2 3 4 5 6 7 8 


a, 

J 

P; 

1 j 


IITTITII 
ITIIITTI 


at 
j 

P7 


3 4 6 

1 3 4 5 8 


at 

P7 


6 4 3 

1 3 4 5 8 




1 3 3 4 4 5 6 8 




3 5 7 


y7 


1 2 4 6 8 



Table 1: Example for the definition of the various sequences needed in the for- 
mulation and in the proof of lemma 1. The spin pattern corresponds to the prod- 
uct D(^i)C(^ 2 )5(^3)5(^4)£ > (^5M(^ 6 )C(^7)£ ) (^8) of monodromy matrix elements, 
m=8, |a+| =3, |P~| =5. 



number may be zero if all monodromy matrix elements . . . 1 are equal to C. We 
claim that the j-fold left action of monodromy matrix elements results in 



M 

<0| Y[C(X k ) T^(X M+l )...T« j (X M+J ) 
L k=\ 



M+Yi M+Y2 M+y„ 



M+j 



I I • • £ cp 1 (£ 1 , Yl )...( Pn (4,Yn)(0| i]Oh). (A.6) 



h=\ h=\ 4=1 



k=l 

k^i\ ,...,4 



where either 



M+y„ 



1 



or 



(pn(4,Y«) = a{X L )c{X M+ln M n ) VI < n > , 



(A.7a) 



1 



(p«(4,Y«) = ^(^J c (^4Am+y„) II /n >r 



(A.7b) 



Here we defined a sequence (y n )™ =1 by arranging all a+ and in non-decreasing 
order (see table 1 for an example). 
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Proof. Let us prove our statement above by induction over j. First of all (A. 6) is 
satisfied for j = 1 by (A.3)-(A.5) and the trivial formula 



M M+l 

<o| Y[c(h) c(W) = (o|]l c (^)- 

L £=l 



M+l 

n 

jk=i 



(A.8) 



Next, let us assume that (A. 6) is true for some j between 1 and m — 1, and let us 
act with T^ J ^ (kM+j+i) from the right. We have to distinguish four possible cases 
(a 7 -+i,P y+ i) + = (2,1), (1,1), (2,2), (1,2). 



The 



case (2, 1). Then Tj^ 1 (X M +j+i) = C(k M +j+i), and 



M 



k=l 



(oi \[\c{x k )\ t;; ( w) . . . t^; (i M+J+l ) = 



M+y l M+y 2 M+y„ M+j+l 

e 1 ••• 1 c P1 (£ 1 , Y1 )...(p n (4,Y„)<oi n c (^)- ( A - 9 ) 

4=1 4=1 e„=\ k=\ 



This is of the same form as (A. 6), and the number of sums n is in accordance with our 
rule, since a, + i =|= 2 and |3 7+ i =f = 1. 

The case (1, 1). Now T^ 1 ^ (X,M+j+i) = A(^M+j+i), the number of sums increases 
by one in accordance with our rule, and y n+ \ = j+l. Using (A. 3), (A. 6) we obtain 



M 



<oi [n c H t s w • • • c ( wo = 

k=l 

M+Y1M+Y2 M+y„ M+j+l 

E E ••• E £ <pi(^i,Yi)---<p«(4,Yn) 

€1=1 4=1 4=1 4+1=1 

h±lx 4^i,.-,4-i4+i^i,-,4 



a (^4+i) c (^+j+l,^4 +1 ) 



M+j+l 



1 



n b {x k M n+l ) 

^4,-,4+i 



M+j+i 

(01 n c (^)= 

k=l 

k^£ u ...,£ n+ i 

M+j+l 



M+Y1M+Y2 M+y n+1 

EE" E (Pi(^i,Yi)...(p«+i(4+i,Yn+i)(0| ]1 C (^) (A- 10) 

4=1 4=1 4+1=1 fc=l 
h¥*i 4+i^i,-,4 ^i,-,4+i 



which is again of the form (A. 6). We substituted y n+ \ = j + 1 and used (A.7a) in the 
last equation. The cases (2, 2) and (1,2), being very similar, are left to the reader. □ 



Next we set j = m in equation (A. 6). Because of the constraint (A. 2), we have m 
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sums on the right hand side and thus, 

M 

(oi [n c (^)j r P T (wo • • • T£(\ M+m ) = 

k=l 

M+ji M+y 2 M+y m M+m 

EE" E q>i(^i,Yi)-"9™(^™,Y™)<0| Il c (^)- ( A - n ) 

i x =\ e 2 =i e m =i k=\ 



The right hand side of this equation becomes explicit after deciding for each func- 
tion q>j if (A. 7a) or (A.7b) applies. In order to do so we note that a unique partition 

({Y + } 5 {Y~}) e P2(%m) exists such that 



a + 



Y Y7 = P7, Yy<YififP7 



a7 



(A.12) 



For an example see again table 1. Now each factor q>j(£j,yj) with jj ^ y^ for j ^ k 
stems either from the application of A and must be replaced by (A.7a) or it stems from 
the application of D and must be replaced by (A.7b). In the former case j e {y + } while 
in the latter j e {y~}. For a pair of equal ys, say y 7 = Jj+i, which stems from the appli- 
cation of an operator B, tyj(£j,yj) must be replaced by (A.7b) and 9/+i(^/+i,Yj+i) is 
determined by (A.7a) (see (A.5)). Hence, j G {y~} and j + 1 e {y + }. Taking the three 
cases together, we see that for every j e {y + } we have to substitute 9j using (A.7a), 
while for every j e {y~} equation (A.7b) applies. Thus, we have derived the following 

Lemma 3. [16] The general left action formula, version 1, 



M 



(o| [n c (^)J r pT (WO • • • t£"(W0 

fc=l 



M+y„, 



M+a 



1 



M+Y! M+Y2 

e e e n«(v c (w+>v n &a . v , 

£ 1 = 1 4=1 £ m =l L/=l V J fc =l 0\hkiM+) 



m-\(V | ' vj j M+m 

n ^ )^ am +P7 ) n Z7T- ^ (oi nw- ( a - 13 ) 

i=1 ■> •/ ; Jk=l "r'-'^/J Jk=i 



M+m 



Let us briefly rest with this formula. Equation (A. 13) is a very compact formula for 
the general left action of a (5 z -conserving) string of monodromy matrix elements on a 
state of the form (A.l). In a sense (A. 13) seems even more generic than the formula 
(39) that we are aiming at, since the coefficient under the sum in (A. 13) is entirely 
defined in terms of the functions b{X,n) and c(A,,//) that determine the 7?-matrix and 
the vacuum expectation values a(k) and d(X) of the (quantum) transfer matrix. On 
the other hand we face a rather complicated dependence of the limits of summation 
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on the matrix indices of the density matrix elements and the explicit appearance of the 
partition ({y + }, {y~}) in the coefficient under the sum. 

Starting from (A. 13) we can now derive (39) by rather elementary manipulations. 
Let us nevertheless sketch the main steps. The notation 



|£ ++ |=card{l,...,M}n{£ Y +}' a = + 1 l , 
\£~ + \ =card{M+l,...,M + m}n{£,+ } K 
|£ + -|=card{l,...,M}n{£ Y7 }7 = - |a+l , 
\£—\= card{M + 1 , . . . ,M + m} n {£ yJ }™Z^ 



(A. 14a) 
(A. 14b) 
(A. 14c) 
(A.14d) 



turns out to be convenient in the calculation. We shall denote the coefficient under the 
sum on the right hand side of (A. 13) as S = SVi,...,4r Then 



M+aj 



Y[a(X e )c(X M+a +,X e ) Y[ 

■7=1 j 



1 



m- or 



Jf+PJ 



n rf (^_)c(x v ,v P7 ) n wrji 



=(-i)" 4 



•m— or 



M 



n M V nn^,>n w 



uv, 



where 



U 



|a+|Y+-l 
7=1 k=l 



■m-|a+|Y,-l 

n ii b (\M) 

j=\ k=\ y > 



(A. 15) 



(A. 16a) 



V = 



\ a I ; 



711 (/. 



n ^yw 7 ) n ^ 



(A. 16b) 



In the second equation (A. 15) we have pulled out several factors and we have intro- 
duced the notation 



a(X) y sh(A, — A,jfc — T[) 



(A. 17) 
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Notice that we did not assume the Xj to satisfy any equation here. In the next step we 
simplify the terms U and V defined in (A. 16). 
Simplification of U. 



U = 



|cc+|j-l yj-j -irm-|a + |;'-l Ij-J 

n n h\ + , h + ) n h\ , \ + ) n n h\ , h- } n , \ + 

j=\k=\ k j k=\ k J JL ;=1 Jt=l r * fc=l 'j '* 



7=1 *=1 



7/ Y * 



•m-|a + |j-l -i rm-|a + ||a + | -i r|a + |y— 1 

n riMH-H-) n nn*(w 



7=1 fc=l 



Y; 



7=1 *=1 



Y^ Yj 



(A. 18) 



Here we used in the first equation that yt — j is the number of y^s smaller than y^. In 
the second equation we combined the second and fourth factor on the left hand side 
into the second factor on the right hand side. Now we have removed the dependence 
on from the upper limits of the products. 
Simplification of V. We first observe that 



m m 



■|0C T | -i rm-|0C T | 

vu u sh(^-^)=sh m (Ti) n^ (+) n v t ] 



7=1 k=\ 



7=1 



;=i 



(A.19) 



with 



at sh(^ + -^-Tl) 



_i j J "V' V «+ ^ 'U m 

' =Tn S V EI — ^— n sn (V + -^), (A.20a) 

7 Sh(^ -^+-Tl) sh(X, Y, 



PJ sh(^_-^ + Tl) 



V =Tn e — — HI — i^T p-r- IT sh(A* _-£*)• (A.20b) 



k+M^l 



Now we have to distinguish different cases. 
(i)£ y+ e{l,...,M} k+M^L+, and 

• ; • i 



v} +) =-nshfa+-b-v) n Mh y+ -w- 



aj-l 



(A.21) 



Jfc=l 



k=aj+l 



(ii) £ y + e {M+ 1, . . . ,M + y Y + = M + at} => fc + M ^ £ y + for all > at, 



and 



at-i 

^ ^n^-^- 1 !) n M^ + -^k). (a.22) 

^ k=\ j k=a++l ' 
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(iii) L- e{l,...,M} k+M^£ y -, and 

'7 0' 



vH= n sh(A, v -^+Ti) n s mv-^) 



(A.23) 



(iv)£ Y - e{M+l,...,M+y Y -=M + pT} fe+M^£ Y -foralU>|3T,and 



(-) 



(A.24) 



*=P7 + i 

Combining (A.19)-(A.24) and using the notation (A. 14) we conclude that 

1 



' m m 

n n sh(^-y 



ocj — 1 



n n sh(v + -^-ri) n -&o 



7=1 



£=aj+l 



l-|CC+|P;-l 

n n sh(x v -^+r,) n ^ 

7=1 *=i J ' *=P7+i ; 



(A.25) 



which is again of a form where the range of the products does not depend on any- 
more. 

Let us now return back to the sum on the right hand side of (A. 13). The coefficient 
under the sum is our 5, equation (A. 15), where U and V are given by (A. 18) and (A.25). 
The summation variables £j range from 1 to M + jj, i.e., 

l/v+l 



£ y+ e{l,...,M + a+}, 7 = 1,.. 
l y - e {1,...,M + PT}, ; = l,...,m-|a 



(A.26a) 
(A.26b) 



Let us now assume that £ + e {M + cd" + 1, . . . ,M + ra} for some j. Then /V/ + = 

'7 ' ' Yj 

Z,e+ £ {^a++i' • • • »^m}. and i,: follows from (A.25) that 5 = 0. Similarly, 5 = if 



£„- G {M + 13 ■ + 1 , . . . , M + m\ for some j. Hence, we may extend the upper limit of 

summation to M + m in all sums on the right hand side of (A. 13), since this just means 
to add zeros to the sum. 

Then the sum is symmetric in the summation variables and may be written as 



M+m 

£\ ,...,£ m — 1 
tj&k for Jfr 



M+m 

(oi 

k=\ 



(A.27) 
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Since the product over Cs is symmetric in the £j as well, we may now permute the 
summation variables in S according to our convenience. Substituting 



Jj ' ~|a+|-y+l> 



L+ - > ( 

tyj — > ^\a+\+j, 
we can simplify U, 



j= l,...,|a + |, 
j = l,...,m— |a + |, 



(A.28a) 
(A.28b) 



U 



■m-\a + \j-\ -i rm-|a + ||cc + | 

n m\ II II"> 



■| a+ l 7-1 

nn^Vi-wVu+i' 

7=1 ^=1 



J] b(ke k ,hj), (A.29) 

l<j<£<m 



and 5 becomes 



5=(-sh(Ti)) 



\e++\+\£+- 



M 



■7=|a+|+l -I 1-7=1 k=\ °\ A k,^£j 



m m 



k^£j 
1 



\<j<k<m J L/=l Jt=l SfcJ 

r |cx+| «+-l 

nn^Vi-m-^-i) n *(Vi-^-w 

7=1 *=1 yt=a+ + l 

m-|a+|P7-l 

n nM\ +l+j -^ k +r\) n sh ^,„ +l+ ,-^) 

■ 7=1 &=1 fc=Pj + l 



Here we introduce 



= P7 



;= l,...,|a + |, 
j= |a + | + l,...,m 



(A.30) 



(A.31a) 
(A.31b) 



and shift the indices of the products in the last two terms in (A.30). Setting \£ + \ = 
\£ ++ \ + \£ + ~\ and substituting S back into (A. 27) we obtain the desired result (39). 
Notice that, because of (A. 14a) and (A. 14c), we have 



|£+|=card{l,...,M}n{£ 7 }7 =1 
in the new summation variables (A. 28). 



(A.32) 
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Appendix B: The coefficients in the generating function 

In this appendix we give a proof of theorem 2. For fixed n e {0, 1 , . . . ,ra} consider the 

nth term £>a!^a™(^i, • • under the sum (71). It satisfies oci H hOC m = 2m — n, 

which means it has |a + | = n up-spins. Since the sequences of upper and lower indices 
of D are identical, it follows that {5+}^ =1 n {^j}J =n+i = 0. Hence, ({6c+}, {{J - }) G 
P2(Z m ), and a unique P e & m exists, such that 



Pat = 7, 7=1,. ..,n, P|3 y =7, 7 = 71+1, 
Let us define a sequence (£/)™ =1 by 

-i ;'=!,.. .,7i 



(B.l) 



+ i 7 = n+l, 



,m. 



(B.2) 



Then 



n n ^(o);-^-Ti) n 

7=1 fc=l £=0C+ + 1 

m P; — 1 m 

n n sh(Q> y --^+Ti) n sn K-^) 

7=n+l fc=l *=P7+1 

= J| sh(o)p k -^j + 2ep k r\)sh(o) P j-^ k ) (B.3) 

l<j<k<m 

which considerably simplifies the integral formula for D"|;;;[^ (£i , . . . , £ m ) . 

We wish to express the product Yll<j<k<m s h((0yt — K) y' + r l) as wei l m terms of (Op/. 
For this purpose we note that the range M of the multiplication variables consists of all 
ordered pairs (j,k) and is, in general, not invariant under P. However, if (£,n) e PM , 
then (ti, ^ P;7tf . It follows that 

Yl sh((o yk -(o ; + ri) = sh(sign(Pk-Pj)(o) Pk -0) P j)+r\) 

l<j<k<m \<j<k<m 



CO, - COjfc 



sh((0^-Q)p J + sign(Pfc-P7')ri) 

l<y'<^<m 

= sign(P) Yl M & Pk - iOpj + 2e Pk r[) . (B.4) 

l<j<£<m 



Here we used that the first product in the second line is a ratio of two van der Monde 
determinants. We further employed the identity 



sign(Pfc - Pj) = 2e Pk , for j < k, 



(B.5) 
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which can be inferred from the definitions of P and (ey). Inserting (B.3) and (B.4) into 
(34) and performing the Trotter limit we obtain 



n ai...a m /£ e 
^tti ...ouASl) • • • ) S« 



f = \Jc 27ii(l + o(to y -)) 



n 



,j=n+l Jc 



2jii(l + a(co/)) 



sign(p) n 

\<j<k<m 



sh((D Pfc - ^ + 2z Pk x\) sh((D Pj - & 
sh((Dp k -0) Pj + 2£ Pk r[) 



det(-G((D ; ,y) 

(B.6) 



This is a nice alternative form of the integral formula for the density matrix elements 
in the special case of only diagonal local operators e\ and e\. 

The reason why we introduced the seemingly more complicated expression (B.4) 
into (B.6) is that in this form the invariance properties with respect to certain varia- 
tion of the permutation P become obvious. Namely, if Q G & m with Q{l,...,n} = 
{!,... ,n}, then 



n ai...a m /£ 
iy ai...a m ^Sl 



ft/ ^ fr / 

) = Vc 2m(l + a(co;))JL,iVc 



c 2jui(l + a(©/)). 



det(-G(a)y,^)) 

Ill</<iKm Sh(^--4k) 



sign(2P) \\ 

\<j<k<m 



sh (&Q Pk - & Q p j + 2£ QPk r\) 



(B.7) 



where we have used that, by definition of the permutation Q, the identity ZQp k = £p k 
holds for all k G {l,...,m}. 

Next we would like to carry out the summation in (71). The sum can be equiva- 
lently written as a sum over all partitions ({6c + }, {|3~}) G p2(JL m ) with |oc + | = n. As 
is clear from (B.l) every such partition uniquely corresponds to a permutation P G & m 
which hence will be denoted P{a+}- It follows that 



> • • • > Sm) — 



2ra(l + a(co y -)) 



.JJ^c 27ii(l + a(co y -)) 



det(-G((o ; ,y) 

Ill<7</t<m s h(^y — ^) 



£ si g n(ep {6t+} ) n 

({Sc+},{p-})e P2 (Z m ) l</<*<m 



(B.8) 



for all Q G & m which leave the set {l,...,n} invariant. These 2 form a subgroup, 

n, then /'{g+j 7^ ^{5+'}, and 



a 



say il c & m . If now {oc+} ^ {a +/ } but |a+ 
QP{a+}{^ + } = {!,• ••,«} but QP{a+>}{<*> + } 7^ {!,- ••,«}, which means that fl 
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HP{a+/}. = 0. Thus, il determines a decomposition of the symmetric group & m into 
cosets, 



6 m = 



u 



HP, 



{a+}- 



(B.9) 



({5+},{p-})ep 2 (Z m ) 
|cc+|=n 



With this information we can average equation (B.8) over the subgroup il which 
has n\(m — n) \ elements, 



/n(^l) • • • 5^m) — 



n 



J(0, 



./=Vc 27Ui(l + a(CQ;)) 
det(-G(co;,^)) 1 



...rijic 27ii(l + o(co y -)) 



ni</<JK™sh(^-^) n\(m-n)\ ^ X si § n (2 P {5+}) 

in<j<k<m vbj w v ) Sell ( { 5+},{p-})ep 2 (Z m ) 

|a + |=n 

n sh ( (0 e/ J {5+} ^ ~ \i + 2e e p {5+ ^ri) sh(a) 6 p {5+}J - - jU) 



.yU^c 2m(l + o(co 7 -)) 



det(-G(a);,^)) 



n 



!<./<fc<m M1 lS/ 



} = Vc 27ti(l + o(co y -)) 

- 1 - £ sign(P) n sh ^-^ + 2£ ^) sh (^-^ . (B.10) 
n ! (m - n) ! p fg m l </<^< m sh (a> Pt - + 2e PJt Tl ) 



It is now easy to establish a relation with the main result of [26]. Define 



A* = ooifc + ejfcTi . 



(B.ll) 



Then 



n!(m - n) ! p £g m i</<l< m sh(tO Pi - Q) Pj - + 2^) 

= 1 y . , p v T-r sh(V fc - + EpfcTj) sh(Vj -% k - Z P fC\) 

n\{m-n)\ p |g m Slgn !</<!<„ sh (Vfe - V 7 + (e P * + e Py )Ti) 



= G„(m|{X}|{^}). 



(B.12) 



The latter equation defines G„(m|{?i}|{^}) which agrees with the function defined in 
equation (3.8) of [26] upon the identification r\ = i£. Using the formulae in theorem 
3.1 and above lemma 3.2 of [26] and inserting back (B.ll) we arrive at the desired 
result (75). 
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